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ABSTRACT 



We combine extinction maps from the Two Micron All Sky Survey (2MASS) with Hipparcos and Tycho parallaxes to 
obtain reliable and high- precision estimates of the distance to the Ophiuchus and Lupus dark complexes. Our analysis, 
based on a rigorous maximum-likelihood approach, shows that the p-Ophiuchi cloud is located at (119 ± 6) pc and the 
Lupus complex is located at (155 ± 8) pc; in addition, we are able to put constraints on the thickness of the clouds and 
on their orientation on the sky (both these effects are not included in the error estimate quoted above). For Ophiuchus, 
we find some evidence that the streamers are closer to us than the core. The method applied in this paper is currently 
limited to nearby molecular clouds, but it will find many natural applications in the GAIA-era, when it will be possible 
to pin down the distance and three-dimensional structure of virtually every molecular cloud in the Galaxy. 

Key words. ISM: clouds; dust, extinction; ISM: individual objects: Ophiuchus complex; ISM: individual objects: Lupus 
complex; Stars: distances; Methods: data analysis 
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1. Introduction 

In recent years, much effort has been dedicated to the study 
of dark molecular clouds and their dense cores. One of the 
main motivations for these investigations is the study of 
the process of star and planet formation in its entirety, and 
a deeper understanding of the effects of the local environ- 
ment. A key aspect of the scientific analysis of a dark molec- 
ular cloud is its distance, which is related to many physi- 
cally relevant properties (in particular, the mass scales as 
the square of the distance) . Unfortunately, the distance es- 
timates for many clouds are often plagued by very large 
uncertainties and it is not rare to see in the literature mea- 
surements that differ by large factors. Often, this quantity 
is inferred by associating the cloud to other astronomical 
objects whose distance is well known. For example, the 
Lupus complex is located near the Sc o 0B2 association 
(|Humphrevslll978l : Ide Zeeuw et aLlfTQQQl ) , whose distance is 
estimated_to_be ^ 150 pc based on the photometry of OB 
stars ()Krauuteij|199lD . However, this method is based on 
some degree of arbitrariness when making the link between 
the cloud and the other objects, and thus the deduced dis- 
tance can be completely unreliable. 

The cloud complexes considered in this paper, the p 
Ophiuchi and the Lupus dark clouds, are good examples of 
how different distance estimates made by different authors 
can be. Quoted values for the p-Ophiuchi clou d are in the 
range 1 20-16 5 pc, a s recently summarized by IRebull et al.l 
(|2004D . [Ctiiiiil (|1981[ ) estimated the distance to be ~ 160 pc 
from multi-color photometry of heavily absorbed stars in 
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the p Oph core, while iKnude fc Hod (|l998t) suggested a 
significantly smaller figure, ~ 120 pc. For Lupus, the sit- 
uation is even more extreme, wit h estimates from 100 p c 
(|Knude fc Hod Il9981 to 190 pc (jWichmann et al.lll998[ ). 
both based on Hipparcos data; the most widely accepted 
distance is in th e middle of these two extremes, at approx - 
imately 140 pc (jHughes et al.lll993l ; [de"Zeeuw et al.lll999D . 



These uncertainties severely hamper our understanding 
of the physical properties of molecular clouds, and thus our 
knowledge of star formation. For example, an error by a fac- 
tor two on the distance of a cloud translates into an error 
by a factor four in the mass, and has a thus a huge impact 
on the estimate of the density, size, stabi lity, and star for- 
mation efficiency of their cloud cores (see lAlves et al.ll2007l 
for further discussions on this point). 



In this paper we present accurate distance measure- 
ments of the Ophiuchus and Lupus complexes based upon 
Hipparcos and Tycho parallaxes. Our technique is statisti- 
cally sound and, when applied to nearby giant molecular 
clouds, is able to provide accurate distance measurements 
and related errors. The paper is organized as follows. In 
Sect. [5] we present the main datasets used and discuss a 
simple approach to the problem. A more quantitative tech- 
nique is developed in Sect. [3l where we also presents the 
results obtained for the two clouds considered in this pa- 
per. The technique is further developed in Sect.|3]to include 
the effects of the cloud thickness and orientation. Finally, 
in Sect. [5] we briefly discuss the results obtained. 



2 M. Lombard! et al: Hipparcos distances of Ophiuchus and Lupus cloud complexes 




Guliiclic Longitude 

Fig. 1. The Ophiuchus region considered. The gray map shows 2MASS/Nicer extinction, with overlaid the Ak — 
0.1 mag (smoothed) contour. The dots indicate the position of the Hipparcos stars used for the analysis. The boxed area 
indicates the central part of the cloud. 



2. Basic analysis 

Our technique is based on a simple fact: stars observed 
through high column-densities must exhibit a significant 
reddening, while foreground stars will show no reddening 
even when observed in projection toward dense regions of 
the cloud. The first step was thus to make reliable extinc- 
tion maps of the p-Ophiuchi and Lupus clouds in order to 
delineate the cloud boundaries and regions to be consid- 
ered for the distance analysis. For the purpose, we used 
approximately 42 million stars from the Two Micron All 
Sky Survey (2MASS) point source catalog to constru ct a 
1 672 square degrees Nicer (iLombardi &: Alved 1200 iD ex- 
tinction map of the Ophiuc hus and Lupus dark neb ulas. The 
map, described elsewhere (jLombardi et al.|[2008D . a has a 
resolution of 3 arcmin and a 3cr detection level of 0.5 vi- 
sual magnitudes. We considered the two regions shown in 
Figs. [T] and [21 approximately corresponding to the galactic 
coordinates 



-12° <l<2° 
for Ophiuchus, and 

353° < I < 345° 



9° < 6 < 24° 



4° < 6 < 19° 



(1) 



(2) 



for Lupus, and select ed there all Hipparcos and Tycho 
(jPerrvman et al.lll997t ) parallax measurements of stars in 



regions characterized by a relatively high column densi- 
ties (more precisely, in regions with Ak > 0.1 mag in 
our 2MASS extinction map). Finally, in order to estimate 
the intrinsic colors of the Hipparcos stars, we obtained 
their spectral types. For the pur pose, we used the All-Sky 
Compiled Catalogue fASCC-2.5: TKharchenkolbOOlD . which 
lists 2 501 304 stars taken from several sources (mainly 
the Hipparcos- Tycho family, the ground-based Position 
and Proper Motion family, and the Carlsberg Meridian 
Catalogs), and the " Tycho-2 Spectra T ype Catalog" 
(Tycho-2spec, HI/231; IWright et aP MM), which cross- 
references 351 863 Tycho stars with spectral catalogues 
(mainly the Michigan catalogs). In summary, our joined 
catalog contains for each star observed in projection with 
relatively large extinction regions the parallax and its re- 
lated error, the spectral type, and a set of visual magnitudes 
with errors. 

We then considered all Hipparcos stars with well de- 
termined parallaxes (we required the parallax error to be 
smaller than 3 mas) in the region considered, and esti- 
mated the intrinsic B ~ V of each star from its spec- 
tral type (for this purp ose, we used the color tables in 
lLandolt-Bornsteinlll982l p. 15 , and the normal reddening 
law of iRieke fc Lebofskvlll985D . Finally, we evaluated the 
extinction of each star by comparing its observed B — V 
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Fig. 2. Same as Fig. [1] but for Lupus. 



color with its estimated intrinsic color. As final selection, 
we required stars to have a y-band extinction error smaller 
than 1 mag, which left 180 objects in the Ophiuchus region, 
and 186 objects in the Lupus region. A plot of the star col- 
umn density versus the Hipparcos parallax for both cloud 
complexes is shown in Figs. [3] and El a distance-reddening 
zoom plot for Ophiuchus is presented in Fig. 2) 



Figure[3]proves the effectiveness of the parallax method: 
from the right to the left, stars are initially found close 
to the Ay = mag line, with a relatively low dispersion; 
then, as we go to smaller parallaxes, an impressive "wall" of 
stars with significant reddening is observed (approximately 
at TT = 8 mas) . Since for Fig. [3] we only selected stars ob- 
served against high-column density regions in Ophiuchus, 
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Fig. 3. The reddening of Hipparcos and Tycho stars ob- 
served in direction of high integrated column densities 
{Ak > 0.1 mag) in the Ophiuchus field. Only stars with 
relative error on the parallax smaller than 0.3 were plot- 
ted here. In this and other plots the reddening has been 
converted into T^-band extinction by assuming a normal 
reddening law (jRieke fc Lebofskvlll985f ). 



S 1 




Distance (pc) 

Fig. 4. The reddening as a function of the Tycho dis- 
tance for Ophiuchus stars. The plot shows the same points 
as Fig. [3l but is done in distance instead of in paral- 
laxes (and, as a result, error bars in distance are approxi- 
mated). Filled dots are stars for which the local extinction 
is Ak > 0.15 mag, open dots stars with Ak > 0.1 mag. 
This figure can be better used to obtain a precise distance 
of the cloud, that we evaluate as d = (120 ± 10) pc. 



and we further excluded stars with suboptimal measure- 
ment errors, the parallax-reddening plot appears extremely 
clear and gives alone an accurate measurements of the dis- 
tance of the cloud, that we estimate to be approximately 
120 pc. 

The same plot for the Lupus region, shown in Fig. [51 
appears to be slightly less clear, for several reasons. First, 
Lupus has a complex structure and is composed of several 



Parallax (mas) 

Fig. 5. The reddening of stars as a function of their par- 
allaxes in the Lupus region (filled circles: Ay > 1.5 mag; 
open circles: Ay > 1 mag). Stars with relative parallax 
error larger than 0.5 were excluded from this plot. If the 
cloud were a perfect "wall" of dust, we would obtain a 
Heaviside function with the discontinuity at the parallax 
of the cloud, much like Fig. [3] for Ophiuchus. Measurement 
errors, the complex structure of the cloud, and the small 
number statistics all make this plot to more scattered and 
thus less clear. 



subclouds spanning several degrees (see Fig. [5]); hence, it is 
not unlikely that different clouds are located at different dis- 
tances (see discussion below). In addition, the Lupus cloud 
complex seems to be at a larger distance than Ophiuchus, 
at the limit of the Hipparcos sensitivity, and hence the plot 
of Fig. [5] is perhaps not as crystal clear as the one of Fig. [31 
Still, a preliminary qualitative analysis indicates a cloud 
distance close to 150 pc for Lupus. 

An intrinsic problem of this technique is that, in this ba- 
sic formulation, it does not allow a simple, robust estimate 
of the error associated with the distance measurement. In 
addition, the value provided must be regarded as an upper 
limit, since the distance is typically estimated from the first 
star that shows a significant reddening. 



3. Likelihood analysis 

A less qualitative estimate can be obtained by using a sta- 
tistical model fo r the p arallax-reddening relation. Following 
iLombardi et al.l ()2006f ). we assumed that the Ay measure- 
ment for foreground stars is small, while a fraction / of 
the background stars is distributed as a normal variable 
with relative large mean and variance. Mathematically, we 
considered the two distributions 



p'^iAv) = Ge.u{Av\A'^,af^ 



p^s^Ay) = Gau(v4y + a_ 



(3) 
(4) 



where we denoted with Gau(a;|a;, cr^) the value at x of a 
normal probability density with mean x and variance a^. 
The two distributions p^^ and p^^ refers to foreground (i.e., 
almost unextincted) and background (heavily extincted) 
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Fig. 6. The likelihood function of Eq. ^ for Lupus (left 
black histogram) and Ophiuchus (right black histogram) as 
a function of the cloud parallax TTdoud marginalized over 
all the other parameters. Both histograms shows marked 
peaks at ttlup — 7 mas and Troph — 9 mas, corresponding to 
distances of approximately 119 pc and 155 pc, respectively. 
The grey histogram refers to the likelihood function for the 
central region of Ophiuchus (marked window in Fig.[T|). 



stars. These are taken to be normal distributions with vari- 
ances given by the sum of the intrinsic scatter of column 
densities (c^^ and and of the measurement error for 

the star considered {a\^). Note, in particular, that cr^^ 
accounts for the local changes in the cloud column den- 
sity with respect to the average value A^, while a^J^ is 
introduced to allow for underestimates of the photometric 
errors on the star magnitudes or for foreground, thin nearby 
clouds. Finally, the probability distribution p( Ay |7r) for the 
Ay measurement for a star at parallax tt was taken to be 



p{Av\t:) 



\p^^{Av) if TT > TTcloud 

[ (1 - /)/s(Ay) + /p'^S(Ay) if ^ < TTeloud 



(5) 



Note that the parameter / G [0, 1] denotes the "filling fac- 
tor," i.e. the probability that a background star is extincted 
by the cloud. 

We considered this simple model leaving all parameters 
{■TTcloud, f,A^.^, a ^^,A\^,a^^} free. In order to assess the 
goodness of a model, we computed the likelihood hmction, 
defined as 



cloud ; 



N 



Abg bg s, _ 



^ ^ /'OO 

rr / 1 7r("))p(^(") I ^(")) dTT , (6) 

where the product is carried over all N stars of Figs. [3] and 



[51 and where p( 



is the probability distribution 
that the n-th star with measured parallax tt'^"^ has a true 
parallax tt^"-'. Note that the integral appearing in Eq. ([6]) 
properly takes into account the possibility that a star with 



a measured parallax tt^"^ < TTdoud (or tt^"^ > TTdoud) is 
incorrectly taken as background (respectively, foreground) . 
For this probability we used a normal distribution: 



(7) 



where air is the estimated error on the parallax. Note that 
the use of normal distributions in Eqs. (ji]), and ([7]) 
makes it possible to perform the integral ^ analytically in 
terms of the error function erf. 

The likelihood function © was analys ed using Monte 
Carlo Markov Chains (MCMC; see, e.g. iTanneii Il99l[ ). 
Figure [S] shows the likelihood as a function of the cloud 
distance, marginalized with respect to all other parameters. 
From this figure we can trivially evaluate the confidence re- 
gions for the estimated distances; in particular, we obtain 
rfoph = (119 ± 6) pc, and dmp = (155 ± 8) pc, both at 
the 68% confidence level (since the likelihood function is 
well approximated by a Gaussian, other confidence bounds 
can easily derived from the errors quoted here). We stress 
that the errors quoted here include only statistical errors, 
as deduced from the parallaxes and reddening uncertain- 
ties; other effects, such as the thickness of the cloud or its 
orientation in the sky are not taken into account and are 
discussed below in Sect. |31 We stress also that the Lupus 
distance estimate approximately refers to the average loca- 
tion of background stars in the field (i.e., to / = 339° and 
b = 7°), approximately 5.5 degrees south than the center of 
the field considered, and that only a few stars are observed 
through Lupus 1. This difference can play a significant role 
in case of noticeable tilts in the cloud orientation or in case 
of multiple components at different distances (see below). 

Note that the reliability of the maximum-likelihood 
technique presented here was also checked with numerical 
simulations. It is well known that the maximum- likelihood 
method is, under certain circumstances (verified in our 
case) , asymptotically unbiased; however, we decided to test 
its behaviour in our specific case, with a limited (although 
large) number of data. For the purpose, we simulated a 
"thin" cloud together with observation of a population of 
stars. We included reasonable errors on both the star par- 
allaxes and column densities, and recovered the distance 
of the cloud using the same method described here. These 
simulations confirmed that our estimate is essentially un- 
biased (i.e., that on average it returns the true distance of 
the cloud), and that the error provided is rehable. 

We also considered a two-screen configuration, with 
two partially overlapping clouds at different distances, and 
checked the distance provided by the method. Interestingly, 
the simulation showed that in this case the method tends to 
return a "weighted" distance, i.e. a weighted mean of dis- 
tance of the two screens (where the weights are essentially 
proportional to the number of stars intercepted by the two 
screens and by the column densities associated with the two 
clouds). This result is very encouraging, since it guarantees 
that, at least in the relatively simple cases considered, the 
maximum likelihood provide a sensible and unbiased esti- 
mate of the cloud distance. Note, however, that the last 
result also shows that in case of composite clouds with 
significantly different densities of background stars (e.g., 
because of their different galactic latitude), the maximum- 
likelihood technique will be biased toward the distance of 
the cloud with the highest density of stars. This point is 
particularly relevant for Lupus, a cloud complex that spans 
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several degrees in galactic latitude, and for which there are 
indications of possible different distances among the vari- 
ous s ubclouds (see below Sect. 14.21 cf. also iAlves &: Francol 
120061) . 



4. Second order geometrical effects 

Given the angular size of the regions considered, it is not 
unlikely that the cloud complexes have a thickness of 15 pc 
or more. Similarly, even for "thin" clouds, we can easily 
imagine that they are not perfectly perpendicular to the line 
of sight and that different regions of them are at different 
distances. All these geometrical effects can be introduced 
in the maximum-likelihood analysis discussed above. 



4.1. Cloud thickness 

The thickness of a molecular cloud can be investigated 
by modifying the probability distribution of Eq. This 
can be done in several ways, which corresponds to differ- 
ent "definitions" of the thickness of a cloud. We chose to 
parametrize the thickness as a continuous, linear transition 
in the filling factor /, which effectively becomes a function 
of the parallax. Hence, we wrote 



p{Av\7t) = /(7r)/s(A^) + [1 _ f{n)]p's{Av) , 
where 



(8) 







if TT > TTl , 



/max(7r - 7ri)/(7r2 - TTl) if Tr2 < TT < TTl . (9) 
, /max if TT < 7r2 . 

In other words, stars up to the lower distance of the cloud 
(tt > TTl ) show no extinction, and thus their column density 
is distributed according to p^^{Av); similarly, stars back- 
ground to the cloud (tt < 7r2) follow, with probability /max 
(the usual filling factor), a distribution p^s(^j{yy^ finally, 
stars embedded in the cloud (7r2 < tt < tti) present an ef- 
fective filling factors that raises linearly within the cloud 
from to /max. Clearly, both parallaxes tti and tt2 were 
taken as free parameters. Note that, similarly to Eq. (O, 
the use of a simple linear function for /(tt) guarantees that 
we can express analytically the integral of Eq. ((6]). 

The results of a MCMC analysis of the likelihood func- 
tion ^ with the distribution of Eq. ([5]) are shown in Fig.[7| 
for Ophiuchus and in Fig. [H] for Lupus. The plots in these 
figures indicate that, unfortunately, the quality and quan- 
tity of the available parallaxes is not sufficient to put strong 
constraints on the thickness of the clouds. In particular, for 
Ophiuchus we obtained a median thickness of 28l^g pc, but 
the errors are clearly very large; the situation is even more 
extreme for Lupus, with a thickness estimate of 51^35 pc. 
Still, this analysis is important because it can in principle 
provide different distributions for the distances of the two 
clouds (see top plots in Figs. [7| and [8|) with respect to the 
ones obtained from the simpler analysis of Sect. [3] (Fig. [6]). 
In this respect, the fact that the two results obtained are 
almost identical confirms the robustness of our distance es- 
timates. In addition, the much larger thickness estimate 
for Lupus and the flat plateau evident in the right plot of 
Fig. [5] suggest that the apparent thickness of Lupus might 
be the result of different Lupus subclouds being at differ- 
ent distances (see below Sect. 14. 2| ). Finally, we note that the 




100 120 
Distance (pc) 

Fig. 7. A density plot of the likelihood of Eq. ([6|) with for 
Ophiuchus with the extinction probability distribution ([8|) , 
as a function of the cloud distance and thickness (the like- 
lihood has been marginalized with respect to all other pa- 
rameters). The three curves enclose regions of 68%, 95%, 
and 99.7% confidence level, while the line-filled corner in 
the top-left indicate a region that is physically excluded 
(observer embedded in the cloud) . The two plots at the top 
and at the right show the distance and thickness marginal- 
ized probabilities, together with the median value (solid 
line) and the usual three confidence regions (dashed and 
dotted lines). 




150 

Distance (pc) 



Fig. 8. Same as Fig. [7] for Lupus. 



method proposed here to study the cloud thickness will cer- 
tainly provide much more interesting results when a large 
statistics of high-quality parallaxes will become available. 



4.2. Cloud orientation (tilt) 

In order to investigate the presence of possible tilts in the 
orientation of molecular clouds, we considered a simple 
model of a thin cloud located at a distance that is a linear 
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Xtilt 

Fig. 9. A density plot of the derived probability distri- 
bution for the tilt parameters dx and dy of Eq. (fTU]) in 
Ophiuchus. Similarly to Fig. [7l the three curves enclose re- 
gions of 68%, 95%, and 99.7% confidence level and the two 
plots at the top and at the right show the X tilt and Y tilt 
marginalized probabilities. 



function of the sky coordinates: 

d{x, y) =do + dxX + dyy . (10) 

In this equation, do is the distance of the cloud at a refer- 
ence point (typically, the center of the region considered), 
and the sky coordinates (x, y) are taken to be measured in 
pc and oriented along the horizontal and vertical axes of 
Figs. [T]and[2] (so that dx and dy are pure numbers). This 
simple prescription can be directly introduced in Eq. ([5]) 
without any further significant modification of the model. 

We studied again this model with MCMC and obtained 
the probability distribution for {dx,dy) shown in Figs. [5] 
and llOl In summary, for Ophiuchus we marginally detected 
a small tilt oriented such that the Ophiuchus streamers 
would be closer to us than the p Ophiuchi core. Indeed, 
if we repeat the simple analysis of Sect. [3] to the window 
marked with a solid line in Fig. [U the denser central regions 
of Ophiuchus that include the well known p Oph core, the 
distance estimate increases to (128 ±8) pc. In contrast, for 
Lupus we detect a vertical gradient in the distance, with 
the high-galactic latitude regions at larger distances than 
the low-galactic latitude ones (note that, according to the 
right plot of Fig. [TOl the gradient in galactic latitude is 
detected at a Str level). Hence, this result together with 
the large extent of Lupus on the sky seems to indicate 
that this complex might not be co mposed of physically 
rela ted clouds, as also prop osed by iBertout et al.l (|1999[ ) 
and lKnude fc: NielsenI (j200l[ ). Unfortunately, in the current 
form the method presented here cannot automatically dis- 
tinguish multiple components of a cloud complex located 
at different distances. Of course, it is always possible to re- 
peat the analysis in different subclouds, but this requires 
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Xtilt 

Fig. 10. Same as Fig.[S]for Lupus. 



an a priori assumption on the subcloud division. For Lupus 
subclouds at high galactic latitudes, in particular, this is 
manifestly not viable given the relatively small number of 
stars present there. As described above, we applied this 
technique to the central region of the Ophiuchus complex, 
obtaining a somewhat higher distance for it. We also con- 
sidered subregions of Ophiuchus occupying different sectors 
of Fig. [T] (defined by the two diagonals, the central verti- 
cal line, and the central horizontal line). However, within 
the relatively large errors of these analyses, we could not 
identify any significant difference in the distance. Note that 
this result is consistent with large errors shown by the plot 
of Fig. [9l which shows that a "flat" screen molecular cloud 
is still consistent with our data. 

4.3. Effects of diffuse ISM 

So far, in our analysis we have implicitly assumed that the 
only significant extinction is associated with the molecular 
cloud. In reality, the dust present in the interstellar medium 
(ISM) is known to produce an average y-band extinction 
of approximately 1 mag kpc""'^ along the galactic plane. In 
principle, it is not difficult to take into account this effect 
by replacing Eqs. ^ and ([4]) with 

p'HAvIt:) = GMAv\A'^+c^/7r,cTX+alJ , (11) 
p'^^Av) = Gau(Ay + a/TT, + J , (12) 

where a ~ 1 mag mas. In practice, however, given the rel- 
atively large scatter tr^^ of Ay compared to a/n for the 
cloud considered in this paper, the effect of the ISM can 
be safely neglected for the foreground star probability dis- 
tribution p^s. In reality, we argue that the same conclusion 
applies also to the background distribution p^^. First, we 
note that the large majority of Hipparcos stars have rela- 
tively large parallaxes: for example, in the Ophiucus region 
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we found only 10 objects out of 180 with parallaxes larger 
than 2 mas. In addition, the combined effect of a large scat- 
ter in cr^^ of and of the filling factor / makes the de- 
tection of a gradient in A^^ currently unrealistic. A specific 
maximum-likelihood analysis using Eqs. (Ilip and (jl2p has 
been carried out, and the results show that the arguments 
just presented hold. In particular, no significant diff'erences 
in the distance determinations were found using either a 
fixed 1 mag mas or a variable value for a. 

We stress, however, that with the advent of GAIA, it 
will be possible to apply the method described in this paper 
to significantly more distant clouds, for which the cumula- 
tive effect of the diffuse ISM can play a relevant role. 

5. Discussion 

Our analysis indicates that the Ophiuchus cloud is very 
likely to be closer than the we think: the "standard" dis- 
tance, 160 pc, is indeed excluded at Scr. Not surprisingly, 
ou r estimate is in excelle nt agreement with the one obtained 
bv lKnude fc Hod (|l998l ) using a similar technique (but note 
that the two results are based on a different datasets and 
use different statistical met hods). 

Recently, IWilsonI (|2002f) and lMamaieS (|2007t ) reported 
an Hipparcos estimated distance of Ophiuchus of d = 
136^7 pc and d = 135ty pc, respectively, results which are 
marginally in a greement w ith our distance estimate of the 
Ophiuchus core. lMamaiekl mainlv focused on the Lynds 168 
cloud, and inferred the cloud distance from the parallaxes 
of seven stars presumably associated with illuminated neb- 
ulae. Although his method is plagued by a potential uncer- 
tainty, the real association between the Hipparcos stars and 
the Ophiuchus cloud complex, and suffers from the small 
number statistics (7 parallaxes), the agreement obtained 
is reassuring and suggests that we finally have in hand a 
reliab le estimate of the Ophiuchus cloud complex. IWilsonI 
(|2002f ). instead, used a more rigorous maximum-likelihood 
approach, but his analysis is substantially different from 
ours. In particular, his method is based on a preliminary 
classification of stars as foreground and background based 
on their observed reddening; the distance of the cloud is 
then inferred from the constraints imposed by the paral- 
laxes of the stars considered. In contrast, in our method we 
do not need to explicitly classify stars, but on the contrary 
can use directly all data in the likelihood (an advantage of 
this is that we are able to use in our method also stars that 
cannot be uniquely classified as foreground or background 
because, e.g., of thei r large photom et ric err ors). 

Very recently, iLoinard et al.l (|2008f ) used phase- 
referenced multi-epoch Very Long Baseline Array (VLBA) 
observations to measure the trigonometric parallax of the 
Ophiucus star-forming region. Their method is based on 
accurate (to better than a tenth of a milli-arcsecond), ab- 
solute measurements of the position of individual magnet- 
ically active young stars though to be associated with the 
star-forming region. If multi-epoch data are available, one 
can infer the parallax of the star, and thus of the star- 
forming regions, with exquisite accuracy (below 1 pc for 
objects within ~ 200 pc). This technique was applied to 
four stars in the p-Ophiuchi region, and interestingly it lead 
to somewhat contradictory results: two objects associated 
with the sub-condensation Oph A (SI and DoAr21) have a 
measured parallax close to 120 pc, and thus in excellent 



agreement with our estimate; two others, associated with 
the condensation Oph B (VSSG14 and VL5), have a much 
higher distance of ~ 165 pc. A possible explanation of this 
is provided by a closer examination of the two sources in the 
Oph B condensation. From their spectral energy distribu- 
tions (SEDs), they appear to be Class III A- type stars, with 
little or no circumstellar emitting dus t, and an extinction 
as large as 50 mag i n the V band (see IWilking et aIlll989L 
and especially Green e et aLlllQQ^ for a discussion on the in- 
dividual SEDs of VSSG14 and VL5). These results indicate 
that VSSG14 and VL5 are probably background stars, not 
directly associated with the Ophiucus complex, and likely 
to be part of a smaller OB association at ^ 165 pc. 

The distance of the Lupus complex is substantially in 
agreement with the ones reported in the literature, but both 
the thickness analysis and the orientation analysis suggest 
noticeable differences in the various Lupus subclouds. 
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